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We propose the experimental realization of (3 + 1) relativistic Dirac fermions using 
ultracold atoms in a rotating optical lattice or, alternatively, in a synthetic magnetic 
field. This approach has the advantage to give mass to the Dirac fermions by coupling 
the ultracold atoms to a Bragg pulse. A dimensional crossover from (3 + 1) to (2 + 1) 
Dirac fermions can be obtained by varying the anisotropy of the lattice. We also 
discuss under which conditions the interatomic potentials give rise to relativistically 
invariant interactions among the Dirac fermions. 
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For their high level of control, trapped ultracold atoms are ideal systems for simulating 
in a tunable way strongly interacting models jlj]. A well-known example is the experimen- 



tal realization of interacting lattice 
transitions has been both detected 



iamiltonians: for bosonic gases, the Mott-superffuid 
2j and investigated in a variety of interesting situa- 



t-ups jl|; f< 

studies 3), 4] of metallic and insulating phases of a two-species mixture in a 3D optical 
lattice have opened the way to experimentally investigating the rich phase diagram of the 
Fermi-Hubbard model. 

The examples mentioned above refer to the ability of cold atom systems to simulate 
non-relativistic Hamiltonians but a fascinating new challenge of the field is the tunable 
experimental realization of relativistic systems which are relevant to high energy physics 



and quantum gauge theories 



the properties of graphene 



5|. It is worth mentioning, for instance, the simulation of 



6(, i.e. (2 + 1) relativistic Dirac fermions, obtained by using 



ultracold fermions in honeycomb lattices 



7 



10( . Other recent proposals to realize massless 



(2 + 1) Dirac fermions consist of ultracold fermions on a square lattice coupled with 
properly chosen Rabi fields 



111, interacting bosons in a two-dimensional lattice produced 

n 

by a bichromatic light-shift potential with an additional effective magnetic field [12[ and 



13|- 



bosons with internal energy levels in a tripod configuration 

It is then highly interesting to see whether it is possible to go beyond the (2+1) case 
and simulate relativistic (3+1) Dirac fermions. We are concerned, in particular, with the 
possibility to make them massive and also interacting, possibly in a Lorentz invariant way. 
Mixtures of two ultracold fermionic species (and recently of three species 

HQ) 

may 

also be useful for the experimental realization of Dirac fermions with internal degrees of 
freedom. New developments in this direction could open the way to simulate, by cold atom 
systems, Kogut-Susskind staggered lattice fermions jltl 1?J or more general elementary 
particle theories. In perspective, this development could permit to study in a controllable 
experimental set-up part of the phase diagram of QCD |18| . 

The aim of this paper is to discuss an experimental scheme to realize (3 + 1) massive 
Dirac fermions (with a mass eventually time-dependent) using ultracold atomic fermions, 
a set-up which makes possible to control interactions through Feshbach resonances [lj] and 
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to realize mixtures of different internal states. A method of simulating the Dirac equation 
in (3 + 1) dimensions for a free spin-1/2 particle in a single trapped ion was presented 



in 



19] . where the transition from massless to massive fermions was also studied. Here 
we propose instead to use non-relativistic polarized ultracold fermions in a rotating cubic 
optical lattice with tight-binding Hamiltonian 



H = -tJ2 ( c ! ( iA , c., + h.c. 



where c] creates an atom in the i-th well of the lattice, t is the tunneling parameter 
(assumed for the moment equal along the three axes x, y and z) and the sum is on 
nearest-neighbours wells. The lattice (created with three counterpropagating laser beams) 
is assumed to be rotating with angular velocity u so that the electrically neutral atoms feel 
an effective magnetic field: with the minimal substitution —iHV — > —iftV — mA (where m 
is the mass of the atoms and A = u x r is the analog of the magnetic vector potential) we 
have the Hamiltonian flTJ, with = [mjK) J! 1 A-dl. Rotating lattices have been efficiently 



realized quite recent 
optical deflectors 



y employing four intersecting laser beams manipulated with acousto- 



21]). Alternatively, one could also end up in the Hamiltonian (P) using, 



d obtained by spatially 



22]. 



on a cubic lattice, fermions subjected to a synthetic magnetic fie 
dependent optical coupling between internal states of the atoms 

Before studying the spectrum of the Hamiltonian (JTJ, let's briefly comment on the 
reason of its choice: one may wonder, in fact, if a simpler Hamiltonian - without a magnetic 
field - of the form H = —tY2ij c l^ij c j (Bij = 1 if i and j are nearest-neighbours and 
otherwise) is able to simulate (3 + 1) Dirac fermions. This is equivalent to ask whether it 
is possible to realize, with a suitable choice of Bij, a semi-metal such that the bands touch 
at isolated points. In 



23] the symmetries groups which lead to a spectrum without Fermi 
surface and energy gap were classified: although this result can be used to exclude certain 
classes of B's matrices, it does not help however to identify the tight-binding Hamiltonians 
which could have the desired spectral properties. By direct inspection, we checked that the 
3D Bravais lattices with a single atom per cell and only nearest-neighbour hoppings does 
not give band touching in isolated points at zero energy. It is for this reason that we focus 
our attention on the realization of Dirac fermions using an artificial uniform magnetic 
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field B = lotA = 7r0o(l, 1, 1), corresponding to u = (h/8ma 2 )(l, 1, 1) (a is the lattice 
spacing and <fto = h/2ma 2 ). With a ~ 2/xm, one gets for K atoms a rotation frequency 
v = \u\/27T ~ 200Hz. From now on we set </>o = 1 and a = 1. 

The magnetic field 5 = n( 1,1,1) induces a 7r-flux on every square face (see the 
schematic plot in Fig.l): to diagonalize Hamiltonian (CEJ) is mostly convenient to use the 
gauge A = 7r(0, x — y,y — x), similarly to 24|. The quasimomenta k take the values in the 



magnetic Brillouin zone [26|], given by —ix/2 < k x>y < n/2, —n < k z < n |27j . 

Since the lattice sites can be divided in two inequivalent sets, say A and B (see Fig.l), 
we can write the Fourier transforms as cAk) = Xljer c j e%k ^ (where T = A,B) and plugging 
in eqn. (CQ), we get the spectrum [24j, [25[: 

E(k) = ±2 1 y/cos 2 k x + cos 2 k y + cos 2 k z . (2) 

An energy spectrum like ([2]) was obtained for PbTe-type narrow-gap semiconductors with 
antiphase boundaries 28[; a model having this spectrum has been recently used in 29| . 
where it was shown that a suitable distortion of tunneling couplings in fermionic lattices 
can introduce a scalar and a Yang-Mills field. 

For half-filling the Fermi energy is zero and there is a vanishing gap between valence 
(E < 0) and conducting bands (E > 0) at the isolated Dirac points k = ±|(±1, ±1, ±1). 
A pair of inequivalent Dirac points is given by kR = |(1, 1,1) and k^ = — 1(1, 1, 1). 
Expanding the energy around these Dirac points we have E^k^m + q)/h ~ wf|<z|, where 



o ii o :: <f V 




FIG. 1: Sets A (red squares) and B (blue circles) on the cubic lattice (right: artificial magnetic 
field B). 
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the Fermi velocity is given by vp = 2ta/ H. Close to these zero-gap points, the quasiparticles 



behave as massless (3 + 1) Dirac fermions of both the chiralities 30J and the linearized 



form of the Hamiltonian (ED) becomes, in the continuum limit, the 3-d Dirac Hamiltonian 

(3) 



H = -lit J df (ip R a ■ Vi) R - rj)\a ■ V^l) 



where a = (cr x ,a y ,a z ) are the Pauli matrices and the two- components spinors ^(f), 
ipR{r) are respectively the Fourier transforms of 

/ (T\ ( c A{k + k R )\ r ( c A (k + k L )\ 

Mk) =\ _ _ ;Mk) = ^ _ • (4) 

\c B (k + k R ) ) \c B {k + k L ) J 

In experimental realizations, the magnetic field B = ir(l, 1, 1) may be subjected to 
some fluctuations which change the magnetic flux per plaquette, 27r$, around the value 
$ = 1/2. In the thermodynamical limit L — > oo, where N = L 3 is the number of sites of the 
cubic lattice, these fluctuations are expected to influence the Dirac cones because, when 
the flux on a plaquette is different from a rational number -, the usual Bloch functions 
are no longer a faithful representation of the translation group and therefore the energy 



spectrum assumes a fractal structure 3lN33|. Note, however, that: (a) experiments with 
trapped ultracold atoms are done with finite number of atoms (N ~ 10 3 — 10 5 ); (b) 
for finite L, the spectrum is not sensibly affected by fluctuations of $ which are much 
smaller than 1/L. To clarify this point, consider two close rational values of the flux, say 
| and j^: in the second case one has q = 100 inequivalent vertices (and q sub-bands), 
while there are only q = 2 in the first case. However, if L 1/5$, the q sub-bands are 
gathered in two groups, each of them almost degenerate. Hence, for the two values of 
the flux we have two physical situations which, for realistic numbers L and 5$ 1, are 
practically indistinguishable. This can be explicitly checked by numerically diagonalizing 
the Hamiltonian ([T]): e.g., even for a rather small size (L = 16) with open boundary 
conditions and a fairly large value of 5$ (~ 5% of $ = 1/2), the spectral density is not 
sensibly affected and it is in reasonable agreement with the one computed from O . 



Let now expose the ultracold atomic gas to a Bragg pulse (see for example 



35|). 
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For deep optical lattices, the Hamiltonian ([T]) acquires a new term of the form 

(5) 



where the sum runs on the lattice sites while k Bragg and u are the differences between 
the wave-vectors and the frequencies of the used lasers. The Bragg term (jSj) gives rise to 
a mass of the Dirac fermions: choosing ksragg = k^ — k R = (tt,tt,7t), the quasiparticles 
around the Dirac point ki (k R ) are transferred close to the Dirac point k R (ki) inverting 
the chirality. This is equivalent to add a mass term to the Dirac Hamiltonian ([3]) 

y cos{ut){if)\if} R + if> R if} L ) = y cos(ut) tpip (6) 

( i)R \ A 1 \ 

where if) — I and if> — if)' I (0 an d 1 are the 2x2 zero and identity 

matrices). Notice that above we exploited the periodicity of the magnetic dual lattice. 
When the frequency difference vanishes, u = 0, one has a time-independent Dirac mass 
while, keeping k Bragg fixed, but changing randomly the intensity of the two lasers, one has 
also the interesting possibility to realize a Dirac fermion with random mass. 

The term ([5]) is actually a particular case of a more general situation: starting from 
the Hamiltonian H = —t ^ ct e~ l ^ Aij+AijS> Ci + J2j Ao(J) c ] c j an< ^ repeating almost un- 
changed the calculation that led to Q, one gets in the continuum limit the Dirac Hamil- 
tonian in an e.m. field 

H = -lit j df (ip f R (r) a ■ D ip R (r) - if} j L (r) a ■ D if) L {rj) 
+ / df A (ip R (r)ip R (r) + ip[{f)ip L (f)J 

where the perturbations Aij and A are slowly varying in space and time and D = V + A. 

So far we have considered equal hopping parameters along the x, y and z— axes but, 
since the tunneling rates depend on the power of the lasers, one can easily realize different 
hopping parameters t x ,t y ,t z . The energy spectrum in this case is 



E(k) = ±2 Jtl cos 2 k x + t 2 y cos 2 k y + t 2 z cos 2 k z . (8) 

and the isolated Dirac points are therefore unaffected by the anisotropy of the hopping 
parameters. The energy spectrum can be also derived including next-nearest-neighbour 
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hopping rates and if they are small with respect to t, as it happens in graphene [6j, the 
low-energy dynamics is still well described by Dirac fermions. 

Using these anisotropic hopping parameters, we can easily realize the crossover from 
(3 + 1) to (2 + 1) Dirac fermions: to this aim, it is sufficient to lowering an hopping 
parameter to zero (say t z , amounting to increase the power of the laser along z) while 
keeping fixed the magnetic field. When t z — > 0, for the fermions we have ipn(k + k^) — > 
<7xipi{k + k' R ) and ip L (k + k L ) -> a x a z ip 2 {k + k' L ), with k' R = (f,-f,-f) and k' L = 
(—§,§>—§), clearly equivalent to q R and q^. The Hamiltonian, on the other hand, becomes 
a purely two-dimensional one 

H 2D = 2tJ dp (tpl D (p) a -p^l D {p) 

where p = (k x , k y ) is the quasimomentum on the x — y plane, a = i(a x , a y ), 70 = a z and 

i/$(P)= lime**' 7 ' ^1,2 (k±(^,±\ 



"K , TX TX 

2' 2 



This is nothing else that the Hamiltonian for (2 + 1) Dirac fermions obtained in 36]. 



ience, in the 2D limit, we obtain directly a pair of (2 + 1) massless fermions, as expected 



37j . One can also show that, in the limit t z — > 0, the Bragg term (jSJ) gives mass also to 



the (2 + 1) Dirac fermions. 

Notice that, in the discussion above, one can consider two (or eventually more) 
fermionic species: they can be either different hyperfine levels of the same fermionic species 
or different species of a mixture (e.g., a Li-K mixture). Experiments with collisionally 
stable mixtures of two jlj] and also three Q, [lsj]) fermionic species has been recently 
reported. The low-energy Hamiltonian will be then simply the sum of free Hamiltonians 
of the type ()3]): in this scheme, the mass of different Dirac fermions obtained by Bragg 
pulses can be in principle different. 

Let now finally consider two-body interactions among ultracold fermions of the form 



eqn describes general non-local interactions among atoms of the same species, as it 
may be realized in p-wave channels 



. We assume for simplicity that Uij is a function 



s 



only oir — i—j. It is well known that the interatomic term (Q does not generally give 
rise to a Lorentz invariant interaction among the Dirac spinors. In our case we point out 
that, to have Lorentz invariance, a necessary condition is 

U(0) = 0, (10) 

where U(k) denotes the Fourier transform of U(r). This condition implies in fact the cor- 
rect inversion of the chirality of the bilinears entering the terms like ip'^ipip. Moreover, 
for the low-energy dynamics considered here, the two-body interaction reduces approxi- 
mately to f/(|A;|) ~ C/(|A|) [where A = k^ — k^ = 7r(l, 1, 1)], resulting in the locality of 
the bilinears in the quartic interaction terms: using the expressions (jlj) in the interacting 
non-relativistic Hamiltonian, one gets 



U(\A\) J d 4 x{^){^). (11) 

Although not strictly renormalizable, the term flTT]) induces an effective interaction term 
(with coupling parameter having dimensions [mass] -2 , similarly to the interacting term 



in the Nambu-Jona Lasinio model 4l|)- The term ( TIT]) is present also when Bragg pulses 
are used to make massive the Dirac fermions. 

Obviously the interatomic term (Q could even spoil the picture in terms of Dirac 



spinors. However, with a weak two-body potential (U 
guarantees the validity of the relativistic description 



r) -C t), the same condition (fTO]) 



391 ] since this holds as far as the 



lifetime r of the quasi-particle excitations is finite and V = h/r smaller than their energy, 



which is indeed our case (the lifetime can be estimated along the lines discussed in 40] 
and the details will be published elsewhere). 

In the 2D limit eqn (fTTl) remains still valid and it induces a coupling between the 
different spinors ipf D and of the (2 + 1) Dirac fermions; moreover, relaxing condition 
ffTQ]) - not strictly necessary in the 2D limit - one obtains all the coupling terms with an 
even number of ip\ D or ip% D . 

To conclude, we observe that when M different fermionic species are considered the 
interaction term ([9]) reads Yli,j 12a p=i ^tf c i;a c i;a c ].p c j;P- 111 general, the interaction 
intra-(a=/3) and inter-(a^/3) species are different: however, the condition (11 01) . neces- 
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sary to have Lorentz invariant interaction terms, simply reads U a,l3 (0) = 0. If the in- 
teraction is independent on the internal degrees of freedom a, ft, then eqn (TIT]) becomes 
U(\A\) J d 4 x (^ a )(^)- 

Conclusions: We have shown that ultracold atoms in a rotating optical lattice are able 
to simulate (3 + 1) Dirac fermions, with their mass generated by a Bragg pulse which 
transfers particles from a Dirac point to the other. When the two lasers of the Bragg 
pulse have the same frequencies, the Dirac mass is time-independent, otherwise one has a 
sinusoidal time-dependence of the mass. This property could be used to study adiabatic 
or quenched dynamics in the Dirac equation; with random Bragg pulses, it can be used 
to investigate instead diffusion and disorder in relativistic quantum mechanics. 

We have also analyzed the crossover from (3 + 1) to (2 + 1) Dirac fermions which can be 
induced by anisotropic lattices. Finally, we have also given a criterion for the interatomic 
interactions in order to get relativistically invariant effective interaction terms. Interesting 
perspectives along this line include the possibilities: (i) to study the relativistic Hamil- 
tonian for several species with general interactions (in particular with no intra-species 
interactions); (ii) to simulate the Nambu-Jona Lasinio model 41] ; (hi) to manipulate the 
ultracold atomic lattice for realizing Majorana fermions; (iv) to study the (eventually 
attractive) relativistic interacting theory which we have obtained, also in the 3D — 2D 
crossover. 
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